SHORT RANGE SHORT-TERM 



FUMIGATION MODEL 



FOR THE INCO SUPERSTACK 



Report # AQM-AMB 1 



NOVEMBER, 1980 



TD 

883.7 

.06 

V46 

1980 




Ministry 
of the 
Environment 



The Honourable 
Harry C. Parrott, D.D.S., 

Minister 

Graham W. S.Scott, Q.C., 
Deputy Minister 



1I0W 



Sudbury Environmental Study 
Coordinator: E.W, Piche 
Ontario Ministry of the Environment 
6th Floor, kO St. Clair W. 
Toronto, Ontario 
Canada M*»V 1P5 



Short Range Short-Term Fumigation 

Model 

For The INCO Superstack 



Author 

A. Venkatram 

Atmospheric Model Development Unit 

Ai r Resources Branch 

Ministry of the Environment 

880 Bay St., hth Floor 

Toronto, Ontario 

Canada M5S 1 Z8 



Report Number A °- M ~ mD 



17274 



Copyright Provisions and Restrictions on Copying: 

This Ontario Ministry of the Environment work is protected by Crown 
copyright (unless otherwise indicated), which is held by the Queen's Printer 
for Ontario. It may be reproduced for non-commercial purposes if credit is 
given and Crown copyright is acknowledged. 

It may not be reproduced, in all or in part, part, for any commercial purpose 
except under a licence from the Queen's Printer for Ontario. 

For information on reproducing Government of Ontario works, please 
contact Service Ontario Publications at coDvright@ontario.ca 



Abstract 
List of Tables 
List of Figures 
Acknowl edgements 



ENTERED JUL 2 5 2001 



TABLE OF CONTENTS 



MOF 

STANDARDS DEVELOPMENT BRANCH 

LIBRARY 



Page 

i 

ii 
i i i 

iv 



1. Introduction 

2. The Model 

3. Model Testing 

k. Summary of Model Equations 

5. Comparison of Concentration Estimates with 
Observations 

6. Applications of the Model 

7. Summary 



1 

2 

8 

13 

15 

17 
19 



References 

Tables 

Figures 



20 



Abstract 

This report presents a new model to predict the 
dispersion of pollutants emitted into the boundary layer associated 
with sunny days. The model is applicable to the 1NC0 superstack 
plume which causes air quality problems, primarily during convective 
conditions, when atmospheric turbulence is generated by surface 
heating due to high incoming solar radiation. Under these 
conditions, conventional Gaussian models based on Pasqui 1 I -Gi f ford 
curves cannot be used to estimate ground level concentrations. 

The model has been validated with field data collected 
in Sudbury and elsewhere. The input parameters of the model depend 
on routinely available meteorological measurements. So, it is 
ideally suited for utilization in a supplementary emission control 
system such as the one presently being operated at the International 
Nickel Company in Sudbury. 
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1. Introduction 

The 381 m superstack operated by the International 
Nickel Company (INCO) in Sudbury is the single largest point source 
of S0 2 in North America. On an average, it emits around 3000 tons 
of S0 2 every day. Due to the large buoyancy of the plume, the 
effective emission height is about 800 m from the ground. During 
most of the year the atmospheric boundary layer does not grow to 
this height and the plume is not brought down to the ground. 
However, during the daytime hours of the summer months (June to 
September) the boundary layer associated with surface heating by 
solar radiation grows to heights around 1000 m. The plume emitted 
into the turbulent boundary layer is then mixed down to the ground 
and high S0 2 concentrations can occur. Thus, in terms of local air 
quality, the emissions from the superstack are a problem primarily 
during hot summer days. 

The atmospheric boundary layer in which turbulence is 
largely maintained by buoyancy originating from the surface heated 
by solar radiation is referred to as convective. The vigorous 
turbulence of the convective boundary layer leads to rapid mixing 
of pollutants emitted into it. It is only recently (Willis and 
Deardorff, 1978; Venkatram, 1 980) that we have started to understand 
dispersion in the convective boundary layer. Conventional 
"Gaussian type" dispersion models cannot predict ground-level 
concentrations under such conditions. 

In this report we describe a dispersion model to estimate 
concentrations associated with plumes emitted into the convective 
boundary layer. This model can be used in a supplementary emission 
control system for INCO. 



We might point out that pollutants from stacks over 100 m 
tall are primarily affected by convective turbulence (Venkatram, |g80) 
This suggests that the model is applicable to most major pollutant 
sources in Ontario. Also, as high concentrations from these sources 
occur during convective conditions, the model is a vital component 
of long-term models. 

2. The Model 

Before describing the model it is useful to review some 
relevant features of the convective boundary layer. The turbulent 

velocities a and 0* in the upper part of the PBL (2 2 0.1 2.) are 

v w 1 

proportional to the convective velocity scale w A given by 

w^ = (gH^./^) 1 3 (1) 

where H is the surface heat flux, T is the average temperature of 
the mixed layer and z. is the mixed layer height. Thus, w A is one 
of the controlling parameters for dispersion of pollutants emitted 
at heights much greater than the thickness of the shear dominated 
surface layer. We recall that this surface layer does not extend 
much higher than z = |l_| where L, the Monin-Obukhov length, is 



L = -T o U jV 3 /gkH o (2) 



In (2), u A is the surface friction velocity and k is the von-Karman 
constant. 

As the most energetic eddies in the mixed layer scale with 

z. , the integral time scale x relevant to dispersion is z./w. . We 

i m r 1 * 

also see that T is roughly the time taken after release for pollutants 
to become well mixed through the depth of the mixed layer. This 
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interpretation of T allows us to appreciate the definition of the 

m 

non-dimensional travel distance X, 



X = w^x/z.u (3) 

In (3), x is the distance from the pollutant source and u is the 
mean mixed layer wind. 

With the preliminaries of the preceding paragraphs we 
are now in a position to describe our dispersion model. We will 
find it convenient in our development to deal initially with the 
crosswind integrated concentration 



C Y = / C(x,y,o)dy 



CO 



With CO we can avoid discussion of plume intermi ttency in the 
crosswind direction. 

Figure 1 shows a schematic of the physical system being 
considered. Based on visual evidence of the behaviour of the looping 
plume we will assume that we can separate the small scale spreading 
of the plume about its centerline and the larger scale up and down 
motion. The vertical meandering of the plume is associated with 
the convective updrafts and downdrafts extending through the depth 
of the PBL. As these downdrafts (or updrafts) have relatively long 
lifetimes (Lamb, 1978) it is reasonable to assume that a plume segment 
emitted into a downdraft will remain in it until it impinges on the 
ground. The figure illustrates such a situation. Note that the 
vertical velocity of the plume can be resolved into an upward acting 
buoyant velocity and an opposing downdraft velocity which eventually 
brings the plume down to the ground. As the emission point is taken 
to be well above the shear dominated surface layer most of the 



downward travel of the plume occurs in a region in which the 
velocity (and potential temperature) is virtually uniform. 

It is reasonable to assume that a nonzero concentration 
(crosswind integrated) is detected at a receptor x whenever the 
plume impinges on the ground at a distance x < x . Then if Pj(x) 
is the probability density of plume impingement, the fraction of the 
time f the plume is seen at x is given by 



x_ 

/ P d (x,v)dx 



f = / P.fx.vjdx (5) 

o 



In (5), v refers to variables such as stack height, buoyancy and 
atmospheric turbulence which determine P .. 

It is always possible to associate an average vertical 
dimension o to the collection of plume segments detected at x . 
Then, for a given wind speed u in the mixed layer we can express the 
average plume segment concentration at x as 

c ~ A (6) 

p uo 

where Q. is the emission rate. Note that Q/u is the mass per unit 
length of the plume segment. 

We assume that the ensemble averaged concentration 
can be written as 



cv u r ,o) - 2£ ( 7 , 



Note that (7) implies that we are modeling the concentration time 

series at x by top-hat profiles whose constant magnitude is given 

by (6). In other words, the concentration at x is either zero or 

C . We point out that Davidson and Halitsky (1958) based a successful 
P 
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dispersion model on this idea. 

In order to be able to use (7) we need formulations for f 
and o. As plume impingement can occur at any downwind distance 
ranging from to », a convenient choice for P ,(x) is the lognormal 
distribution. The precise form of Pj(x) is not expected to be 
critical in computing C ^ (x ,o) as f depends on the integral of the 
distribution. Then, we can express f as 

x. 

Jinx - Unm 2 

d(£nx) (8) 



■ f Jinx - Unm 

=1 / exp H 2£n*s ^ 

iTr£ns J L 9 



In (8), m and s are the geometric mean and standard deviation of 

g g 

the lognormal distribution. These parameters can be estimated by 
examining the results of numerical experiments by Lamb (1978). His 
study shows that the plume centerline descends when particles are 
released from an elevated source. The rate of descent of the center- 
line (locus of the maximum crosswind integrated concentration) is 
approximately 0.5 w A . Then, for particles released at a height h , 
the plume centerline impinges at a distance wuh /w A for the source. 
Briggs (1975) had a similar picture in mind when he proposed the 
"touchdown" model for plume rise. According to him the buoyant 
plume rises relative to a region of convective downdrafts. As the 
buoyant plume rise rate descreases with time the downdrafts eventually 
bring the plume to the ground at a distance x. given by the solution 
of the following equation 

aF i/3 x .2/3 w - 

! --i-L - -h (9) 

u u s 

where F is the buoyancy parameter and w is the downdraft velocity 
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which is about 0.5 w w . The first term on the left is the buoyant 
rise and the second term represents the effect of convectlve down- 
drafts. Some preliminary studies (Venkatram 1 980) indicate that the 
constant ot in (9) is * 1.0. 

In view of the preceding discussion we feel that it is 
reasonable to assume that the mean impingement distance m is 
proportional to x. . The model testing, described in a later section, 
indicated that the simplest relationship m = x. was adequate. We 
can estimate s by using Lamb's (personal communication) results on 
the statistics of convective velocities in the PBL. He found that 
the magnitude of the downdraft velocity one standard deviation away 
from the mean (= 0.5 w^) was approximately 0.75 w A (= 1.5 w.) . The 
impingement distance corresponding to 0.75 w^ would be a lower limit 
on possible plume impingement distances. We can then approximate 

s by 
9 

x. (w d =0.5 w,J 
5 g = x!(w d = 0.75 wj ° 0) 

It is interesting to note that for the cases considered, s computed 
from (10) was around 2, the value used in an earlier version of this 
model (Venkatram, 1980). 

For the behaviour of a we propose 



— - X ; small X (I la) 

i 



— « constant ; large X (lib) 

i 

Equat ion (l 1 b) reflects the observation that at large X the vertical 
spread is limited by the capping inversion at z,. 



A plausible interpolation for a between the limits for 
small and large X can be written as 



rri 



1 - exp(-1.5 X) (12) 



where A is a constant to be determined from the subsequent analysis. 
The constant 1.5 in (12) was chosen on the basis of model testing 
described in a later section. Also, it is seen that for X > 3, o is 
within \% of its maximum value. This result agrees with the 
observations of Willis and Oeardorff (1976). 

We can now write the expression for the center line ground- 
level concentration as 

dx.o.o) - « 03) 

uo a 

y 

We can determine the constant B in (13) by noting that the 
expression for C(X,o,o) should reduce to that corresponding to the 
well-mixed PBL for large X 

C(X,o,o) - — & (14) 

/2tF uo z. 

y i 

From (H) we see that B/A = 1//2tt. Note that (\k) implies that the 
concentration distribution in the crosswind direction is Gaussian, a 
description which is adequate according to the tank experiments of 
Will is and Deardorff (1976). 

In a recent study Nieuwstadt (1980) showed that under 
convective conditions a ~ X for X as large as k. Moore (1975) also 

y 

found this linear behaviour in his measurement of a values due to 

y 

elevated sources. Based on these studies we selected the following 
s imple form for a 



a /z. = 0.i»5 X (15) 

The constant multiplying X in (15) was suggested by Willis and 
Deardorff (1978). 

It is useful to note that the expressions for the a's do 
not include the effects of self- induced spread due to plume buoyancy. 
The model testing described in a later section indicated that these 
effects were minor. 

It should be apparent from the preceding discussion that 
each stage of model development relied heavily on comparison with 
concentration observations. Thus, this model is essentially a 
concise description of the observations made under a variety of 
meteorological and stack conditions. We point out that a model is of 
little value if it has to be "tuned" each time the dispersion conditions 
change. As will be seen in the next section, our model has been 
applied to four independent sets of data with no-site specific 
"tuning". 

3. Model Testing 

The model was tested with four sets of data. Two of them, 
reported by Weil (1977), refer to measurements made in the vicinity of 
the Dickerson and Morgantown power plants situated in Maryland. The 
remaining data consists of measurements made around the INCO nickel 
smelter in Sudbury, Ontario. Simultaneous meteorological data were 
also obtained during these field studies. 

For the sake of completeness, we will describe some of the 
relevant information given in Weil's report. The coal-fired Dickerson 
plant consists of three 1 85 MW generated units with two 122 m stacks, 
60 m apart. It is situated in the rolling terrain of Montgomery 



County. Measurements at the Dickerson plant extended from October 
1572 through April 1973. 

The Morgantown power plant is also situated in relatively 
flat terrain next to the Potomac River in southern Charles County. 
It has two 575 HW units exhausting through two 213 m stacks, 76 m 
apart. Measurements at the plant were made during February 1975 
through June 1 975 - 

The INCO (international Nickel Company) smelter is 
situated on flat ground in Sudbury, Ontario. There are a number of 
small lakes in the vicinity of the plant. Waste gases from the 
smelter are exhausted through a 381 m stack. The high sulphur content 
of the local nickel ore results in the emission of approximately 
3000 tons of S0„ per day. The first field experiment to study the 
concentrations associated with this emission was conducted during a 
two-week period in June 1978. A second field study was undertaken 
during the week of August 16, 1979. 

The method of measurement of S0„ was similar in all four 
field studies. An instrumented mobile van was used to measure ground- 
level S0 2 concentrations. Using available roads, repeated passes were 
made across the direction of plume travel. The traverses lasted from 
i to 1 hour and resulted in about 6 to 10 instantaneous concentration 
profiles. Some of the traverses in the Sudbury field study were 
done with an instrumented twin-engined Aztec aircraft flying at tree- 
top level (~ 30 m) . 

Concentration profiles obtained from the traverses were 
combined to form an Eulerian averaged profile corresponding to an 
averaging time of approximately 1 hour. The maximum of this composite 
profile was taken to represent the ground-level centerline concentration. 
These centerline concentrations were compared with predictions from 



our model. Following Deardorff and Willis (1975) we used data 
corresponding to the range 1 ,5 w A < u < 6 w^ to ensure the applica- 
bility of convective scaling. 

As the surface heat flux was not measured, we assumed that 
it was 0.3 times the incoming solar radiation. We believe that this 
is not critical as the heat flux is used in the computation of the 
convective velocity w ft which depends on the one-third power of the 
heat flux. 

Concentration observations were compared with model 
predictions using the following three steps: 

1. Compute percentage of predictions within a factor-of-two 
of the observations. The factor-of-two criterion is 
based on experience with dispersion modeling and is the 
time-honored method of testing model performance, (A 
Position Paper of the AMS 1977 Committee on Atmospheric 
Turbulence and Diffusion Model). Predictions meeting 
this accuracy criterion were further tested using curve 
fitting procedures. 

2. Regress observations on predictions using the following 
equat ion : 

C(observed) = a + bC(predicted) (16) 
The value of the coefficient of determination r 2 is an 
indication of the performance of the model. Furthermore, 
the regression coefficients 'a' and 'b' should not be 
significantly different from their expected values of 
0.0 and 1.0 respectively, 
3- Regress observations on predictions using the equation: 
C(( 



(observed) = a C(predicted) (17) 



- 11 



Equation (17) would be relevant for model testing if 
observations are lognormally distributed about the 
predicted ensemble mean. As observations indicate that 
this is indeed the case (Csanady, 1973) we feel that 
the r 2 associated with (17) is a good test of model 
accuracy. Note that the expected value of 'a' and 'b l 
is unity. 
Table 1 summarizes the results of model testing on three 
sets of data. It is seen that model predictions compare well with 
observed concentrations. More than 80% of the predictions are within 
a factor-of-two of the measurements. The value of r 2 for the 
Dickerson data set is above 0.70. This indicates a high degree of 
association between model predictions and observations. For the 
Morgantown data there is a tendency for the model to underpredict 
slightly although the r 2 is as high as 0.80 for the logarithmic 
correlation analysis. 

For the 1978 Sudbury data, the model does not perform as 
well as it does for Dickerson and Morgantown. This could be related 
to the fact that our assumption of constant downdraft velocity might 
be questionable when the emission point is in the upper part of the 
convective PBL. Recall that the Sudbury stack is 381 m tall. 

Figures 2, 3 and k illustrate the details of the comparison 
between model predictions and observations. In Fig. (2a) we have 
plotted the ratio of C(observed) to C(predicted) against X/X. . It is 
seen that the model performs as well for X/X. > 1 as it does for 
X/X. < 1. This gives us some confidence in our choice of the log- 
normal distribution to "smear" the concentration distribution around 
the mean impingement distance x. . The scattergram of Fig. (2b) shows 
in more detail the excellent agreement between predicted and 
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observed concentrations. 

It is seen from the Table that the averaging times for 
most of the observations are less than one hour. Furthermore, the 
concentrations represent averages derived from about six instantaneous 
values corresponding to ground level traverses. In view of this 
(See Venkatram, 1979 for a detailed discussion) substantial deviations 
of observed concentrations from predicted ensemble means are not 
unexpected. However, as seen in Fig. (4a) , this scatter Is well 
within the acceptable norms (see AMS position paper on this subject, 
Bulletin of the American Meteorological Society). In fact, 80? of 
the data is within a factor of 1.5 of the predictions. 

In performing a regression analysis between observations 
and predictions we have chosen only those data points which meet the 
factor-of-two criterion. This filtering is justifiable, in our 
opinion, only if the rejected data points are a small percentage of 
the observations. Furthermore, this process ensures (approximately), 
that the expected variances of the observed concentrations about the 
predictions are comparable in magnitude. Note that linear regression 
is valid only if this condition is satisfied (Neville and Kennedy, 
196M. Figure (3a) indicates, as seen before, that the model 
slightly underpredi cts concentrations for the Morgantown data set. 
However, better than 65% of the predictions are within a factor of 
1.5 of the measured concentrations. The high degree of association 
between model results and observations is also evident in Fig. (3b). 

The details of the meteorological and source conditions for 
the Dickerson and Morgantown concentration observations are given 
elsewhere (Venkatram, 1980). In this paper we only present the 
relevant information for the observations made in Sudbury, Ontario. 
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Note that (Table 2) u < 6 w, for all the data indicating that the 
PBL is dominated by buoyant production of turbulence. This also 
justifies the use of convective scaling (Deardorff and Willis, 1975) 
In the analysis of dispersion. 

Figure kb illustrates the comparison between observed and 
predicted concentrations. Considering the expected variances in 
measurements the comparison is gratifying. In addition to the r 2 
being as high as 0.6, the best fit values of 'a' and l b' are very 
close to their expected values. 

Table 2 gives the relevant information for the 1979 field 
study in Sudbury. Note that the mixed layer winds are small compared 
with those measured during the 1978 experiment. This, as we shall 
see later, has implications with regard to our ability to predict 
observed concentrations. An examination of the Table Indicates that 
9h% of the model estimated values are within a factor-of-two of the 
observed values. However, it is seen that the r 2 value is only 0.1. 
Thus it would seem that although the model meets operational require- 
ments, it does not satisfy statistical tests. However, we wi 1 1 see 
in the next section that linear correlation analysis (r 2 ) is not 
always appropriate and in many cases cannot tell whether or not a 
model is working. 



4. Summary of Model Equations 

At this point it is useful to summarize the equations on 
which our model is based. 

C(X,o,o) = 1— - (18a) 

/2tF o o* 

y 

C = Cz. 2 u/Q; X e w A x/z.u (18b) 



]k - 



a /z. = O.hS X 
y ' 



■[■ 



a/z. =|1 - exp(-1.5 X) 



f(X) = 






- / exp(-t 2 /2)dt 
£ J 

-00 

£n(X/X.)/J!.ns 



X. 



w^x./z.u 



(18c) 
(I8d) 
(I8e) 

(I8f) 
(I8g) 



In (I8f), the mean impingement distance x. is the solution of 



pl/3 2/3 _ + h u = 
i d i s 

w, =0.5 Wj. 



x. (w = 0.5 w A ) 
i d " 

x.(w d = 0.75 w A 



(19a) 
(19b) 

(19c) 



The results from an earlier paper (Venkatram, 1980) indicate that 
using a constant value s = 2.0 yields acceptable model predictions 

It is useful to point out that the expression for the 
ground-ievel concentration when the PBL is well-mixed is simply 



C(x,o,o) = 



0-9 0. 
w A z.x 



(20) 



Note that (20) is independent of the wind speed u. In an earlier 
paper (Venkatram, I98O) we showed that this simple expression can 
provide acceptable estimates of concentrations even for distances as 
small as X = 0.*t. However, the acceptability of a concentration 
prediction depends on the particular application of the model. For 
example, estimates better than a factor-of-two might be necessary for 
supplementary emission control during "adverse" meteorological 
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conditions. However, as we shall see in the next section, there is 
a limit to our ability to predict observed concentrations even if we 
have a "perfect" model. 

5. Comparison of Concentration Estimates with Observations 

Even if we can make reliable estimates of concentrations, 
the nature of the convective boundary layer presents us with problems 
when we attempt to compare these estimates with measured concentrations, 
We have to realize that computed concentrations correspond to 
ensemble averages while measured concentrations represent averages 
over a limited set of "concentration" events. It is possible to show 
(Venkatram, 1979) that the expected deviation of the observed average 
concentration from the predicted ensemble mean can be written as 

2(r - i)t f 

e 2 ■ j 1 (21) 

In (21), e is the deviation as a fraction of the ensemble mean 
concentration, F is the ratio of the instantaneous peak to the ensemble 
mean at the receptor under consideration, x F is the Eulerian time 
scale governing dispersion and T is the averaging time (sampling 
time to be precise). In the convective boundary layer T F can be 
approximated by z./u so that (13) reduces to 

2(r - i)z. 

e 2 * = (22) 

uT 

Close to the stack V can be as large as 100 (see Gifford, 1959). The 
ratio increases away from the centerline to relatively large values 
at the edge of the plume. At large downwind distances V approaches 
an asymptotic value of 2. Eq. (22) implies that comparison of 
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observations with predictions in the vicinity of x. is not straight- 
forward. If we take T = 10 at x. (Pasquill, 197*0, u = 5 ms and 
z. = 1500 m we find that £ = 1.22 for a typical averaging time of 
one hour. In view of this, it might be necessary to average over 
several observations before any meaningful model validation can be 
done. It should be pointed out that this problem with comparison 
with observations is related to the relatively large horizontal 
length scales {~ z.) found in the convective boundary layer. In other 
words, the large convective eddies and wind speed determine the 
duration of a concentration event recorded at a receptor. Thus, for 
a fixed averaging time, the number of events N p which constitute the 
observed average can be approximated by 

n e . T f (a) 

To get an idea of a typical value for hL let us take I = 1500 m and 
u to 5 ms . From (23) we find that the observed concentration 
averaged over one hour (T = 3600 s) corresponds to only 12 concen- 
trations events. This tells us that we are likelv to have problems 
in predicting short-term concentrations during convective conditions. 

It is noted from (22) that the deviation between model 
predictions and observations increases as the wind speed decreases. 
This could explain the relatively poor r 2 for the 1979 field study in 
Sudbury. To understand this better, consider the expression for r z 



2 _ i unexplained variance , ? , , 

variance of observations 



If we denote the average of the observations by C, the unexplained 
variance of a "good" model would be approximately e 2 ^. Then, even 
if the model is working r 2 will be low if the variance of the 
observations is comparable to e 2 ^ 2 . Note that e 2 is determined by 
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the inherent nature of turbulence and has nothing to do with 
measurement or modeling errors. It is readily seen from the data 
presented in Table 3 that the expected variance e 2 ^ 2 between obser- 
vations and model estimates is comparable to the variance of the 
observed concentrations. Thus, the low value of r 2 is expected 
from theory. Our analysts indicates that r 2 is a valid indicator of 
model performance only if the expected variance between observations 
and predictions is much smaller than the variance of the measured 
concentrat ions. 

We see from (22) that e 2 is inversely proportional to the 
averaging time T. This means that model performance will improve 
with averaging time. However, we should note that the magnitude of 
T is limited by the time scale of the evolution of the planetary 
boundary layer. Beyond T ™ 1 hour nonstationari ty effects will 
contaminate observations. To illustrate the role of averaging time 
in determining model predictability we have plotted the ratio 
C(observed)/C (predicted) for the 1979 Sudbury data against time of 
traversing which corresponds roughly to T. In spite of the scatter 
in the data, there is a clear trend indicating improvement of model 
predictions with averaging time. 

Our discussion clearly emphasizes the problems related to 
predicting short-term concentrations in the convective boundary 
layer. The implications of this will be discussed in more detail in 
the next section. 

6. Applications of the Model 

The obvious application of the model is In supplementary 
emission control which works as follows: 
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Meteorological inputs to the model are predicted a few 
hours ahead of time. Then the model is used to determine the emission 
rate which will ensure that the maximum expected concentration 
(predicted) will be below the local air quality standard. Cutbacks 
in emission are implemented if an exceedence is likely. The success 
of this type of emission control depends upon a) the ability to 
forecast meteorological variables and b) the validity of the model 
used to predict the concentration. We will confine our discussion to 
po i n t ( b ) . 

A validated physically correct model Is necessary for 
effective emission control. The convective dispersion model 
described in this report is clearly suitable for this application. 
However, as we pointed out earlier, a model prediciton will deviate 
from the observed value. Furthermore, this deviation is large under 
low wind speed conditions when the maximum concentrations are likely 
to be at their highest. This suggests that a model can provide 
guidance rather than the prediction required for emission control. 
Specifically the model tells us that the most important meteorological 
variables are z. , w A and u. The model prediction based on these 
variables will provide information on the range of values in which 
the observed concentrations will lie. The precise use of this 
concentration range will depend on the economics of the SMpplementary 
emission control program. What we want to emphasize is that even if 
we could forecast meteorological variables accurately, we cannot have 
a perfect emission reduction program. There are bound to be a 
certain number of exceedences of the air quality standard. 

Another application of the model is in the prediction of 
long-term averages. In this case, there is good reason to believe 
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that deviations between observations and predictions will be 
minimized in the averaging involved in computing the long-term (annual) 
average. This filtering of the errors will occur only if the short- 
term sub-model in the long-term model is capable of predicting the 
correct ensemble mean. The model presented in this report satisfies 
thi s cr i terion. 

7 . Summary 

This report presents a model to describe dispersion of 
pollutants emitted from a point source into a convective boundary 
layer. The model is applicable to the INCO plume which causes high 
ground-level concentrations primarily during convective conditions on 
sunny days. We should point out that conventional Gaussian models 
based on the Pasqui 1 1 -Gi f ford curves do not work under these conditions . 
The model described here represents a new approach to dispersion 
model ling. 

The model uses some of the most recent understanding of 
dispersion in the convective boundary layer. Furthermore, model 
predictions have been compared with field data collected in Sudbury 
and elsewhere. The results of the comparison are very encouraging. 
The model is simple to use as it depends on routinely available 
meteorological measurements. The mixed layer height z. can be 
estimated from morning temperature soundings and surface temperature 
forecasts. It can also be predicted with relatively simple models 
such as the one described by Carson and Smith (WO. Similarly, the 
convective velocity scale w A can be readily related to the solar heat 
flux or the mixed layer height using formulations described In 
Venkatram (1978). The model is suitable for operational use in a 
supplementary emission control program such as the one used by INCO. 
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TABLE 1. Summary of model testing for the Dickerson, Morgantown and Sudbury data sets. Numbers in 
parenthesis are the standard deviations. 



Data 
Set 

Name 



No. of 
Data 



Points 



Percentage 

Meeting 

Factor-of- 

two 



criterion 



Mean Value 
of C(obs)/C(pred) 
(f i 1 tered data set) 



Linear 
a(ugm ) 



Logar i thmic 



Dickerson 


26 


85 


0.98(0. 3*0 


0.70 


-0.15 


0.98 


Morgantown 


22 


86 


1.20(0.37) 


0.65 


18.86 


0.87 


Sudbury 


25 


80 


1.02(0.27) 


0.53 


9.95 


1.00 



0.7^ 0.91 1.00 
0.80 2.21 0.8^ 
0.60 0.89 1.20 



tsj 



TABLE 2. List 


ing of 


meteorological and 


stack conditions 


for the 


Sudbury 


concentration measurements 








during June 


1978. Mod 


lei predicts the 


impingem* 


tnt distance x. ; 


and C(pre 


d icted) . 










Run 


DA MO YR 


START 


END 


(mV3) 


, Q u 


X 


z i 


w * 


U -1 


x i 


C(0BS) 


C(PRED) 


C(0BS) 




No. 




TIME 


TIME 


(gs _1 ) 


(m) 


(m) 


(ms-1) 


(ms ) 


(m) 


(ygm"3) 


(ygm~3) 


C(PRED) 




1 


10 06 78 


13:45 


14:32 


2082 


28213 


3100 


1040 


2.13 


11.3 


9372 


208 


241 


0.86 




2 


10 06 78 


14:02 


15:00 


2149 


32224 


39000 


940 


2.16 


11.7 


9479 


260 


361 


0.72 




3 


10 06 78 


15:56 


17:12 


2072 


22958 


32000 


740 


2.31 


10.2 


7650 


494 


373 


1.32 




4 


11 06 78 


16:10 


17:11 


1897 


22958 


3000 


1540 


2.80 


13-8 


6987 


182 


241 


0.75 




5 


13 06 78 


13 = 34 


15:00 


2590 


33330 


9200 


1640 


2.58 


10.5 


6994 


520 


593 


0.88 




6 


13 06 78 


13:01 


13:35 


2703 


33884 


15500 


1640 


2.48 


8.9 


6744 


390 


444 


0.88 




7 


13 06 78 


16:30 


17:04 


2481 


27383 


5100 


1840 


2.66 


10.3 


6493 


780 


519 


1 .50 




8 


13 06 78 


16:44 


17:20 


2481 


27383 


18000 


1840 


2.66 


11.5 


7029 


832 


272 


3.05 




9 


13 06 78 


17:23 


17:54 


2489 


30426 


13800 
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2.64 


10.3 


6576 


312 


367 


0.85 


1 


10 
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10:28 


12:43 
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1.93 


9.0 


10543 


572 


552 


1.04 


ro 


11 


14 06 78 


13:03 


13:26 
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2.28 


7.6 
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1.50 


\^j 


12 
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13:31 


14:05 
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1740 


2.35 
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390 


363 


1.07 


1 


13 
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2.54 


5.8 
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0.46 
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5.54 
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0.66 
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0.82 
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1.03 
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16:18 
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2.38 
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1.18 
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1965 
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1240 
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1.22 
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22 
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2327 
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1 .31 




23 
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9.0 
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0.42 




24 
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1.32 




25 


27 06 78 


17:08 


17:39 


2474 


28628 
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1740 


2.58 


8.1 


5762 


416 
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2.33 





TABLE 3. Summary of model testing for field data collected during June 1979 field study in Sudbury. 



Time of 
Sample 

16/08/79 

1336-1 A48 
1458-1525 

1542-1645 



Average 
Time 
(min) 



72 
27 
63 



No. of 
Passes 



9 
4 
8 



x(m) 



5000 
5500 
5500 



a 

(gs->) 



12800 
15300 
20830 



(m) 



1100 

1530 
1300 



w * 
(ms" 1 ) 



2.05 
2.24 

1.93 



u 
(ms-1) 



3-4 
4.1 

4.5 



F C(0bs) 
(mV3) ( U g m -3) 



1920 
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C(Pred) 
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C(Pred) 
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1425-1457 
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1541-1600 
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32 
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67 
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10 
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3 
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1.87 


3-7 


1990 


390 


289 


1.35 
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Fiqure 1. Schematic of physical 
system under study 
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Figure 2a. Variation of the ratio C(observed)/C(predicted) 
as a function of X/Xj. Horizontal lines 
correspond to the selected values of the ratio. 
Data set refers to the Dickerson Power Plant. 
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Figure 2b. 



Scattergram of predicted concentrations against 
observed concentrations. Data set refers to 
measurements around the Oickerson Power Plant. 
Only data meeting the factor-of-two criterion 
have been used. The best fit line corresponds 
to regression of log trans formed concentrations. 
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Figure 3a. Same as 2a except for Morgantown Power Plant. 
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Figure 3b. Same as 2b except for Morqantown Power Plant. 
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Figure 4a. Same as 2a except for Sudbury. 
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Figure 4b. Same as 2b except for Sudbury 
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Figure 5. Variation of C(Obs)/C(Pred) with sampling time 
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